Metal-insulator and magnetic collapse transitions in FeO at high pressure: 

predictions of the LDA+U method 
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The accuracy of LDA + U total energies is tested by predicting the high pressure behavior of 
rocksalt-structured FeO, including metal-insulator and magnetic transitions, equation of state, and 
lattice strain, as functions of the Coulomb repulsion U. For U — 4.6 eV, 6.0 eV and 8.0 eV, two 
self-consistent solutions, one with rhombohedral and one with monoclinic electronic symmetry, are 
found, and for [7=4.6 eV the monoclinic solution becomes more stable than the rhombohedral 
solution at 100 GPa, leading to an insulator-metal transition. With increasing U, metallization 
occurs at higher pressures. Despite the insulator-metal transition, the high-spin magnetic state 
persists to pressures greater than 300 GPa for U > 0. The LDA + U method also gives improved 
agreement with experiments for ground state properties as compared to LDA and GGA methods. 

PACS numbers: 64.70.Kb, 62.50+p, 71.70.Be 
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FeO is a prototypical transition metal oxide (TMO) 
and Mott insulator 0], and the high-pressure behav- 
ior of FeO will help discriminate among possible theo- 
retical models for correlated magnetic insulators. Non- 
stoichometric FeO (wustite) is also important in geo- 
physics as an end-member of the magnesiowiistite solid 
solution, a major phase in the Earth's lower mantle. At 
ambient conditions, FeO has a rhombohedrally-strained 
Bl structure, which transforms to the cubic rocksalt 
structure above the Neel temperature, 198 K M, where 
the spins disorder. At high temperatures and pressures 
(P>120 GPa and T>1000 K), FeO transforms from the 
rhombohedrally-strained Bl phase to a NiAs (b8)-anti- 
NiAs (iB8) superlattice p, 0, pi r|. At room temperature, 
the cubic Bl phase transforms to the strained Bl struc- 
ture at ~ 15 GPa J7|. The Bl structure is maintained up 
to at least 140 GPa J7|, ^| at room temperature. We con- 
centrate here on the behavior of FeO in the strained Bl 
structure to generate predictions that can be compared 
with experiments on the compression of a Mott insulator 
with this prototypical structure. Recent experiments are 
also in conflict regarding the high pressure magnetic be- 
havior of strained Bl FeO. Mossbauer spectroscopy has 
shown a high- to low-spin transition at 90 GPa M, but 
x-ray spectroscopy indicates that the high-spin state is 
preserved at least to 140 GPa ||. Although at low pres- 
sure FeO contains ~ 4-8% Fe vacancies, we consider only 
stoichiometric FeO, since non-stoichiometry is less im- 
portant in FeO at high pressure |10{ . 

The LDA and GGA p) predict an antiferromagnetic, 
metallic ground state for FeO, but it is actually an in- 
sulator with an optical band gap of 2.4 eV |I3]. This 



arises from the underestimation of the local Coulomb re- 
pulsion in LDA. The LDA + U method E3] is designed to 
remedy this by applying an orbital-dependent potential 
between correlated states, and then subtracting a double- 
counting correction. LDA + U has been well tested at 
zero pressure, and gives insulating behavior in reason- 
able agreement with experiment |14fl . Dynamical mean 
field theory (DMFT) also includes an energy-dependent 
self-energy, but is not yet tractable for the large number 
of computations necessary to understand the behavior of 
FeO under a wide range of compressions, strains, and 
symmetries. Indications are that DMFT should be even 
more reliable [ [15[ , though slower, than LDA + U. 

We use the rotationally invariant LDA + U formal- 
ism [n6[ within the full-potential, linear muffin-tin or- 
bital (FP-LMTO) method, and the GGA for the lo- 
cal exchange-correlation functional. In agreement with 
previous LDA and GGA calculations carried out with 
the linearized augmented plane wave (LAPW) method 
0. H Hil: the LMTO calculation gives a high- to low- 
spin transition at 124 GPa with AV = 9.7% for the cubic 
lattice with the GGA (U = 0). When the rhombohedral 
strain is optimized, however, GGA gives a continuous loss 
of Fe magnetic moment from 80 GPa to ~ 250 GPa, also 
in agreement with LAPW. 

U may be determined from constrained density- 
functional computations p0[ , the most careful study giv- 
ing U — 4.6 eV (2l). We also consider U= 6.0 eV and 
8.0 eV since many estimates are higher than that found 
by Pickett et al. Non-overlapping muffin-tin sphere radii 
were fixed at 1.8 bohr for Fe and 1.5 bohr for O. The value 
of J, 0.89 eV, is also determined by the constrained LDA 
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and is close to the atomic value, 
1.23 eV; small changes in J in this range do not effect 
results significantly. 

We find two self-consistent insulating solutions for both 
unstrained (cubic) and strained lattices. The first is ob- 
tained by constraining the electronic symmetry to have 
the same rhombohedral symmetry as the lattice, and the 
second by reducing the imposed electronic symmetry to 
monoclinic and then perturbing the d-orbital occupation 
matrix. In the rhombohedral case, the spin-down ai g or- 
bital has one electron; the other four spin-down orbitals 
(two e g pairs) are empty and therefore raised in energy, 
forming an insulating gap. The monoclinic solution was 
found by starting with the d-orbital occupancies from the 
GGA calculation on cubic FeO. In this solution, the e g 
state (lower in energy than ai ff in the cubic lattice) splits 
into two singly-degenerate states (e 9 — >a ff +b 9 ), and the 
ai 9 state becomes a s . The e g -derived states are lowest in 
energy, and therefore one of these (a g ) is stabilized, while 
the companion (b g ) state is destabilized along with the 
other a 9 (formerly ai s ) orbital. 

For the unstrained lattice at zero pressure, U = 4.6 eV 
produces a band gap of 0.48 eV for the rhombohedral so- 
lution and 0.13 eV for the monoclinic solution. Strained 
FeO has smaller gaps (0.23 eV and 0.08 eV respectively, 
reflecting the increase in band width with rhombohedral 
strain. We find that the band gap at zero pressure is be- 
tween Fe 3d and Fe 4s states, in agreement with optical 
spectroscopy, which shows weak absorptions beginning at 
about 0.5 eV, assigned to Fe 3d/0 2p^Fe 4s transitions, 
and much stronger absorptions beginning at 2.4 eV, as- 
signed to Fe 3d/0 2p- 
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>Fe 3d transitions |2S 
calculations show somewhat less Fe 3d-0 2p mixing at the 
top of the valence band as compared to the spectroscopic 
results. For the upper and lower Fe 3d (Mott-Hubbard) 
splitting in the unstrained lattice, we find that U = 4.6 
eV gives 1.9 eV for the rhombohedral solution and 1.2 eV 
for the monoclinic solution, while for unstrained FeO, we 
obtain gaps of 1.5 eV and 1.1 eV, respectively. The char- 
acter of the band gap is constant for all values of U, and 
increases approximately 0.25 eV for each 1 eV increase in 
U. The Fe 3d-Fe 3d splitting increases by approximately 
0.8 eV for each 1 eV increase in U, in accordance with 
the model proposed by Zaanen, Sawatsky and Allen Q . 
With U = 6 eV, we obtain a Fe 3d-Fe 3d gap of 2.40 eV, 
in agreement with optical experiments. 

At low pressures, the calculated ground state of FeO 
has rhombohedral symmetry and a rhombohedral strain 
that increases with pressure for all values of U, in agree- 
ment with experiment, as shown in Fig. Pfl. Wc find a 
weak first-order transition from rhombohedral to mon- 
oclinic symmetry at approximately 110 GPa when U = 
4.6 eV, with a very small volume change (AV = 0.3%). 
The minority-spin states for the monoclinic solutions are 
much broader (~0.6 Ry) than for the rhombohedral so- 
lutions (^0.3 Ry), and increased hybridization results in 
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FIG. 1: Energy difference as a function of vofume for LDA + 
U solutions for FeO, relative to the unstrained rhombohedral 
solution, (a) U = 4.6 eV, (b) U = 6.0 eV, (c) U = 8.0 eV. Solid 
line, strained rhombohedral; dashed line, strained monoclinic; 
dash-dot line, unstrained monoclinic. 



a favoring of the monoclinic solution over the rhombohe- 
dral at high pressure. The existence of two self-consistent 
solutions with strain-dependent energies explains the re- 
sults of Mazin and Anisimov J25|, who found the low- 
symmetry state for a highly strained lattice along with 
the high-symmetry state for the unstrained lattice, and 
considered this an anomaly of the LDA + U method. 

In the ground state, rhombohedral solution for strained 
FeO, the band gap closes at approximately 180 GPa with 
U = 4.6 eV, a compression of 38% (Fig. ||). The strained 
monoclinic solution, however, becomes stable at 110 GPa 
(Fig. la; the band gap has already closed in this solution 
at around 60 GPa, so that the insulator-metal transition 
is predicted to occur simultaneously with the symmetry 
change). For all cases, we find that the band gap initially 
increases before decreasing with pressure in the expected 
manner, similar to what was seen in NiO[|14fl. Since the 
band gap is Fe 3d-Fe 4s in character, increasing pressure 
causes the Fe 4s states (unaffected by the application of 
U) to be raised in energy, leading to an initial increase in 
the band gap with pressure. At higher pressures, the Fe 
4s states rise in energy above the minority-spin 3d states; 
band broadening then becomes dominant, and the band 
gap begins to narrow. Higher values of U cause the gap 
to close at higher pressure, and increase the pressure of 
the rhombohedral-monoclinic transition. 

The effect of pressure on U is not known. U is thought 
to be rather pressure independent, being an atomic prop- 
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FIG. 2: Band gaps as a function of pressure for different LDA 
+ U solutions: (a) U = 4.6 eV, (b) U = 6.0 eV, (c) U = 
8.0 eV. Squares, unstrained rhombohedral; circles, strained 
rhombohedral; down triangles, unstrained monoclinic; down 
triangles, strained monoclinic. 



erty, but screening of U might increase with pressure. If 
so, that would reduce the transition pressures. Further- 
more, if screening increases significantly in the metallic 
state, a dramatic drop in U at metallization could help 
to precipitate a magnetic collapse transition coincident 
with or close to the insulator-metal transition. Both An- 
derson model simulations [P6| and dynamical mean field 
theory [£7], |2^] show the growth of a quasiparticle peak 
at the Fermi level with increasing pressure prior to the 
overlap of spin-up and spin-down states, suggesting that 
U decreases with pressure but does not become negligi- 
ble. We do not find find a high- to low-spin transition 
with U > 4.6 eV up to 250 GPa. 

Without spin-orbit coupling, our moments at zero 
pressure (3.41-3.49 fi B for U = 4.6-8.0 eV) are slightly 
higher than experiment (3.32^) p9[ , as expected with 
LDA + U. Previous work has suggested that spin-orbit 
coupling is also important in FeO 1 30 , 3L] , but an orbital 
moment at Fe 2+ would increase the discrepancy with ex- 
periment. In the absence of spin-orbit coupling and for 
U = 0, the rhombohedral strain in FeO splits the spin- 



down t2 g states with ai ff below e g , and the effect of U 
does not induce an orbital moment. If we begin with 
the higher-energy LDA+U solution with e g below ai 9 
and then turn on the spin-orbit coupling, we do obtain 
a converged solution with an orbital moment of 0.56 [Ib 
using U =4.6 eV, essentially in agreement with the re- 
sult of Ref.pO], which gave an orbital moment of 0.88 /is 
for U = 8.0 eV. However, we find this solution to be ~ 
10 mRy higher in energy than the solution without an 
orbital moment, so our computations do not support a 
large orbital moment in FeO. 

We find that the LDA + U treatment gives more accu- 
rate structural properties than does the GGA; the equa- 



tion of state and strain predicted by LDA+U are in 
good agreement with experiment (Table 1). Also very 
encouraging is the excellent agreement with experiment 
obtained for the rhombohedral strain and its pressure de- 
pendence using U = 4.6 eV (Fig. ||). Our ground state 
structure with optimized strain gives a — 58.6°; the value 
extrapolated from experimental strain as a function of va- 
cancy concentration suggests a — 58.4° for stoichiomet- 
ric FeO 0. Larger values of U also result in improved 
agreement with experiment as compared to GGA, but 
not nearly as satisfactory as the results obtained with 
U = 4.6 eV. This suggests that the observed rhombohe- 
dral strain and its value with increasing pressure are also 
strongly dependent upon electron correlation. 

The present results represent the first comprehensive 
study of the applicability of the LDA+U method at pres- 
sures above 100 GPa. An upper bound of 110 GPa is de- 
rived for the proposed transition from a rhombohedrally 
strained Bl structure to a monoclinic structure at high 
pressure, which is coincident with the transition from an 
insulating to a metallic state. No high- to low-spin transi- 
tion in the magnetic moment at high pressure is predicted 
with this model, although a strong dependence of U on 
pressure or metallicity could precipitate such a transi- 
tion. Conductivity experiments on FeO at high pressure 
should help to clarify the picture presented here, partic- 
ularly with respect to the variation of the band gap with 
pressure. Since many models can give insulating states 
for FeO at zero pressure, high pressure behavior will be 
a good testing ground for this and other theories of tran- 
sition metal oxides. 
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TABLE I: Equation of state parameters for the universal 
(Vinet) equation of state |32] for the high-spin states of FeO 
obtained with the GGA and the LDA+U methods for cu- 
bic and rhombohedrally strained lattcies compared with ex- 
periment ||. For LDA+U, U = 4.6 eV, J = 0.89 eV. R, 
rhombohedral symmetry; M, monoclinic symmetry. Volumes 
correspond to one formula unit of FeO. A constant 5390 Ry 
was added to total energies. 
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FIG. 3: Optimized strain for FeO as a function of volume. 
Diamonds, GGA/LAPW; left triangles, GGA/LMTO; circles, 
U = 4.6 eV; up triangles, U = 6.0 eV; down triangles, U — 
8.0 eV; squares, experiment. 



of FcO available prior to publication. 
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